Anisotropic Ordering in a Two-Temperature Lattice Gas 
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We consider a two-dimensional lattice gas model with repulsive nearest- and next-nearest-neighbor 
interactions that evolves in time according to anisotropic Kawasaki dynamics. The hopping of 
particles along the principal directions is governed by heat baths at different temperatures. The 
stationary states of this non-equilibrium model are studied using a generalized dynamical mean-field 
medskipgj^olyglg In the stable state the particles form parallel chains which are oriented along the direction 
of the higher temperature. The resulting phase diagram is confirmed by Monte Carlo simulations. 

PACS numbers: 05.50.-fq, 05.70.Ln, 64.60. Cn 
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Systems driven into a non-equilibrium steady state 
can reveal a number of features not characteristic to the 
equilibrium state. An example is the well-known two- 
dimensional lattice gas model with attractive nearest- 
neighbor interaction whose degenerate ground states vi- 
olate the X ~ y symmetry of the system. Below a criti- 
cal temperature the disordered (high temperature) phase 
segregates into a high- and a low-density phase and the 
strip region of the condensed particles can be oriented ei- 
ther horizontally or vertically if periodic boundary con- 
ditions are imposed. There are several different ways 
to provide non-equilibrium conditions and to violate the 
X — y symmetry. In driven lattice gases introduced by 
Katz et al. the particle jumps are under the influence 
of a uniform external electric field whose direction is cho- 
sen to be parallel to one of the principal axes . Another 
way to define anisotropic dynamics is the so called two- 
temperature models are investigated extensively . In 
these former models the particle jumps are coupled to two 
different heat baths, that is, the Kawasaki (exchange) dy- 
namics is characterized by and Ty temperatures along 
the X and y axes. In contrary to the previous dynam- 
ics, here there is no a uniform particle current through 
the system. At the same time, an energy transport ex- 
ists from one of the heat baths to the other one driving 
the system into a non-equilibrium steady state. This ap- 
proach can describe the effect of an alternating field (or 
polarized light) on the ordering process [||. 

In the attractive interacting system the symmetry 
breaking is represented by the orientation of the inter- 
face separating the high- and a low-density phases. Sig- 
nificantly different symmetry breaking can be observed 
in the model discussed by Sadiq and Binder which 
is characterized by repulsive nearest- and next-nearest- 
neighbor interactions with equal strength. In the half- 
filled system the particles form alternately occupied and 
empty columns (or rows). That is, here the symmetry 
breaking appears directly (as a bulk feature) in the four- 
fold degenerate ground states. The ordering process in 
this model is well understood and the effect of the 
uniform driving field on the ordering process was also in- 
vestigated It was found that the phases consisting 
of chains perpendicular to the external field become un- 



stable and the stationary state is oriented along the field 
direction. In other words, the chain orientation of the 
ordered phase can be controlled by the application of an 
electric field with suitable direction. Similar effect is sus- 
pected when using an alternating field. This conjecture 
was the motivation of the present work. Now we will 
study the effect of anisotropic jump rate on the order- 
ing process within the formalism of the two-temperature 
model. 

First we introduce the model and investigate by adopt- 
ing a simple dynamical mean-field theory [ p^ . Linear 
stability analysis has been performed to distinguish the 
stable and unstable solutions. In the light of these results 
the phase diagram is evaluated at higher level (4-point 
approximation) of dynamical mean-field theory. Finally 
the results are compared with Monte Carlo (MC) data. 

We consider a half-filled, lattice-gas model with 
isotropic nearest- and next-nearest-neighbor repulsive in- 
teractions with equal strength on a square lattice. As 
usual the coupling constants as well as the Boltzmann 
constant are chosen to be unity. The kinetics is governed 
by Kawasaki dynamics characterized by single particle 
jumps to one of the empty nearest neighbor sites [ pd] |. 
To avoid the difficulties arising from the non-analytic fea- 
ture of Metropolis rate, Kawasaki hopping rate is used 
for directions a = x oi y: 



?a(Ai/) = 
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1-f exp(Ai//T„) 
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where A_ff is the energy difference between the final and 
initial configurations. This jump rate is anisotropic and 
satisfies detailed balance at temperature Ty [T^) in the 
vertical (horizontal) direction. 

The equilibrium model [T^ — Ty) undergoes a continu- 
ous phase transition at a critical temperature Tc — 0.525. 
In the ground states the columns (or rows) are alternately 
occupied and empty providing fourfold degeneracy. Fol- 
lowing the notation of Sadiq and Binder the lattice is di- 
vided into four interpenetrating sublattices and the 2x1 
or 1 X 2 long-range order is characterized by the corre- 
sponding average sublattice occupations (z = 1, . . . , 4) 
0. In the half-filled system the 2x1 states consisting 
of vertical chains are given as ni = rt4 = (1 + mx)/'2, and 
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n2 = ns = — mx)/2, where —1 < < 1 is the order 
parameter. Notice that irix can be either positive or nega- 
tive. Its absolute value as a function of temperature may 
be estimated by using the traditional mean-field approx- 
imation. Similar expressions are obtained for the 1x2 
states of horizontal chains, namely ni = n2 = {l + my)/2 
and 713 — ni — (1 — my)/ 2. It is emphasized that the 
four possible states are equivalent {\mx\ = \'my\) in the 
equilibrium system. Evidently, in the high temperature 
disordered phase ttIq = 0. 

In the following simple dynamical mean-field approx- 
imation | [Tc| ] we will suppose that the long-range order 
still exits under the present non-equilibrium condition. 
It means that the anisotropic jump rate results in only a 
difference between the solutions characterized by and 
ruy. 

According to Eq. (|l]) the time evolution of ni is given 
by summarizing the contributions of the jumps modifying 
the particle density in sublattice 1 as 



^ = -2ni[(l - n2)gx{E2i) + (1 - ni)gy{Eii)] 



(2) 



where Eij is the average value of AT? when a particle 
jumps from sublattice i to j. Similar equations can be 
derived for the rest of sublattice occupations. Assuming a 
2 X 1 or 1 X 2 structures (vertical or horizontal chains) the 
stationary solution of Eq. (||) obeys the following general 
form: 
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exp(3r7iQ/2rQ) 
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where a — x ov y respectively. Notice that the jumps 
along the chains do not modify the particle densities 
within the chains therefore the value of the correspond- 
ing order parameter depends only on the perpendicular 
jump rate (or perpendicular temperature). The order pa- 
rameter as a function of perpendicular temperature can 
be determined numerically. This simple approximation 
predicts continuous phase transition with a critical tem- 
perature ri*^^^ = 3/4. 

Obviously, we have two different solutions if ^ Ty. 
In this situation the linear stability analysis is used to 
determine which solution remains stable. This method 
proved to be very useful when investigating the effect of 
the uniform driving field on these ordered structure in the 
same lattice-gas model . Now we restrict ourselves to 
the survey of results without dealing with mathematical 
details given in a previous work | [ic| ] . It is found that the 

ordered structure appears when imn[Tx,Ty) < Tc^^^^ 
and the chains in the stable phases are parallel to direc- 
tion belonging to the higher temperature. As mentioned 
above, the value of the corresponding order parameter 
depends only on the perpendicular temperature. As a 
consequence the chain structure remains stable even in 
the limit case when the parallel temperature goes to infin- 
ity. The resultant phase diagram is illustrated in Fig. ^ 




FIG. 1. Phase diagram predicted by the simplest 
mean-field theory. The arrows indicate the stable chain di- 
rections in the ordered regions. 

The above results imply the possibility of reorientation 
in the low temperature region. Namely, exchanging the 
values of Tx and Ty a particle rearrangement is expected 
as it is observed by MC simulations under the effect of 
a uniform field [|lO|. It is rather difficult to realize the 
two-temperature model, but as we mensioned earlier an 
alternate electric field can cause similar behavior. De- 
tailed analysis of the reorientation phenomenon goes be- 
yond the scope of the present work. Instead, we focus 
our attention to have a more accurate phase diagram. 

The prediction of the above simple mean-field theory 
can be improved by using the generalized mean-field the- 
ory at the levels of 2- and 4-point approximations. In 
these calculations we have to determine the probability 
of all the possible configurations on 2- and 4-point clus- 
ters by evaluating numerically the stationary solution of 
a set of equations which describe the time variation of 
configuration probabilities. Details of such a calcula- 
tion are given in previous papers ||l^ , [l3|| . The 2-point 
(pair) approximation fails to reproduce the phase tran- 
sition even in the equilibrium model {Ty — Tx). More 
precisely, this approximation predicts first-order transi- 
tion. Similar discrepancies have already been observed 
for other sytems p2| , p^ . This shortage of the present 
technique may be corrected by calculating with larger 
cluster(s). 

Exploiting symmetries we have (only) four indepen- 
dent parameters to describe the probability of all the 
possible particle configurations on a square at the level 
of 4-point approximation. In equilibrium this method 
suggests a continuous order-disorder transition and the 
predicted critical temperature is Ti"*^' = 0.566 which is 
close to the exact value mentioned above. 

In the non-equilibrium case {Tx ^ Ty) the prediction 
of the 4-point approximation agrees qualitatively with 
those obtained above (see Fig. ^). However, now the 
critical temperature depends on both temperature. The 
numerical results are summarized in a revised phase di- 
agram as plotted in Fig. ||. A weak anomaly has been 
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found when and Ty are close to Tc ^ . Namely, the or- 
dering temperatures increase first when moving away the 
equilibrium. After a low peak the transition temperature 
decreases toward a constant value, Tc'^^\oo) = 0.94Ti'*^'' 
if the parallel temperature goes to infinity. Hence, in 
agreement with the earlier approximation, the ordered 
state still exists even at infinite parallel temperature if 
the perpendicular temperature is low enough. 



perature (as well as the temperature dependence of order 
parameter) depends slightly on the parallel temperature. 

The possibility of orientation (or reorientation) process 
may appear in other two-dimensional systems belonging 
to the family of Ashkin- Teller models [Iw. The present 
two-temperature approach gives some insight into this 
phenomena even for higher dimensions. 



0.5 







disordered 











0.5 1 1.5 2 2.5 3 
T 

X 

FIG. 2. Comparison of phase diagrams predicted by 
4-point approximation (solid line) and MC simulations 
(squares). 

To check the above predictions a series of Monte Carlo 
simulations p5| ] has been carried out varying the tem- 
peratures Tx and Ty. Using periodic boundary conditions 
the simulations are performed on a square system for dif- 
ferent sizes L = 20, 30, 40, 80, and 200. The size depen- 
dence of data has indicated no serious finite size effects. 
We have determined the time averages of the energy and 
order parameter varying the perpendicular temperature. 
These simulations have justified that the ordering process 
remains continuous for all the parallel temperatures we 
studied. From these data we could determine the value 
of critical temperatures. The results confirm the above 
theoretical predictions as demonstrated in Fig. ||. This 
figure vizualizes that the transition temperature depends 
sligthly on the parallel temperature, in a very similar way 
as suggested by the 4-point approximation and its value 
tends to 0.98ri*^'^^ when or Ty goes to infinity. 

In summary, we have studied the non-equilibrium or- 
dering process in a two-dimensional lattice gas coupled to 
two thermal baths at different temperatures. In the equi- 
librium system the particles form parallel chains oriented 
either horizontally or vertically. Using a simple mean- 
field theory including linear stability analysis we have 
found that the ordered structure appears in the nonequi- 
librium case if the lower temperature is less than a critical 
value and the chains are oriented along the higher tem- 
perature. This picture has been confirmed by dynamical 
4-point approximation and MC simulations. Despite the 
prediction of the simple mean-field theory the sophisti- 
cated methods have indicated that the transition tem- 
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